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Abstract 

The one-dimensional chain of coupled oscillators with long-range power-law in- 
teraction is considered. The equation of motion in the infrared limit are mapped 
onto the continuum equation with the Riesz fractional derivative of order a, when 
< a < 2. The evolution of soliton-like and breather- like structures are obtained 
numerically and compared for both types of simulations: using the chain of oscilla- 
tors and using the continuous medium equation with the fractional derivative. 
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1 Introduction 

Nowadays it become clear that anomalous dynamics and kinetics may be not so pathological as 
it was believed earlier. More and more situations has been reported in the literature supporting 
this ^121 El El El- The immediate causes for these anomalies could be attributed to fractal or 
multi- fractal character of phase space, Levy flights, dynamical traps, or long-range correlations 
which are present in many interesting applications ^ |3| |4| |5j . The desolated once fractional 
order calculus revives again in order to discribe such problems. Indeed, equations which in- 
volve derivatives or integrals of non-integer order appeared to be very successful in describing 
anomalous kinetics and transport [Illfil[7ll8ll9l [Tnirrn[T2] . 

However, fractional equations could be rarely derived explicitly from the equations of mo- 
tion or from a Hamiltonian of the model. More often fractional equations for dynamics or 
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kinetics appear as some phenomenological models. Recently, the method to obtain fractional 
analogues of equations of motion was considered for the bunch of problems related to sets of 
coupled particles or other objects that interact via a power-like potential |131 114j . Examples 
of such systems are one-dimensional chain of interacting oscillators, spins, or waves that can 
be considered as a benchmark for numerous applications in physics, chemistry, biology, etc. 
jl6l 1171 1181 1191 12U1 1351 142j . Transfer from the set of Hamiltonian equations to the continuous 
media equation with fractional derivatives is an approximate procedure. Its formal realization 
can be perfomed using the so-called "transform operator" jl4j . Different applications of the op- 
erator have already been used to derive fractional sine-Gordon and fractional Hilbert equation 
[T!^ . to study synchronization of coupled oscillators and for fractional Ginzburg-Landau 
equation jT^. 

Localized in space structures, such as solitons or breathers are the subjects of a great inter- 
est since a long time. These structures have been widely studied in discrete systems on lattices 
with different types of interactions as well as in their continuous analogues (for a review, see 
[341 \'67\ ). Solitons in a one-dimensional lattice with the long-range Lennard- Jones- type interac- 
tion were considered in |SHI • Kinks in the Frenkel-Kontorova model with long-range interparticle 
interactions were studied in [321 • The properties of time periodic spatially localized solutions 
(breathers) on discrete chains in the presence of algebraically decaying interactions were consid- 
ered in and recently in |26j. Energy and decay properties of discrete breathers in systems 
with long-range interactions have also been studied in the framework of the Klein-Gordon |37j . 
and discrete nonlinear Schrodinger equations [in|. A remarkable property of the dynamics de- 
scribed by the equation with fractional space derivatives is that the soliton or breather type 
solutions have power-like tails. Similar features were observed in the prototype lattice models 
with power-like long-range interactions [351 136 . 42, |31IS]- As it was shown in |141 115j . analysis 
of the equations with fractional derivatives (FD) can provide fairly quickly results for the space 
asymptotics of their solutions. Replacement of the system of particles with power long-range 
interactions by an equation with FD will be called fractional equation approximation (FEA). 
As it was shown in |141 I15j . the FEA appears in the infrared limit when the wave number 
k ^ and some other conditions that will be derived later. At the same time, in many applied 
problems equations with FD can appear as an intrinsic feature of the system (see examples in 
reviews (21 El)- Integration of equations with FD needs specific algorithms which are fairly time 
consuming with not yet well defined errors ,46^. 

The goal of this paper is to study a duality of Hamiltonian dynamics of system of particles 
with power-like interactions with the solutions of equations with FD using the transform opera- 
tor. As a model to be studied we consider a chain of interacting oscillators that can be described 
by the sine-Gordon equation in the continuous limit. 

To incorporate long-range interparticle interactions several nonlocal generalizations of the 
sine-Gordon equation were considered in Refs. |42,.42,.44j using the integro-differential equations 
of motion. The generalization of the standard sine-Gordon equation with space FD was proposed 
and numerically studied in |41j . The form of the fractional Sine-Gordon equation was uu — 
^D^u + sinu = 0, where ^I?" is the Riesz space fractional derivative, 1 < a < 2. This equation 
was proposed as an interpolation between the classical sine-Gordon equation and a nonlocal sine- 
Gordon equation. Here, we consider the one-dimentional lattice of coupled non-linear oscillators. 
We focus on the dynamics of soliton-like and breather-like patterns. Following jl3j . we make the 
transform to the infrared limit and derive the fractional sine-Gordon equation which describes 
the dynamics of the lattice. Using both, Hamiltonian lattice dynamics and the fractional sine- 
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Gordon (FSD) equation, we can compare solutions and demonstrate their similarity as well as 
to discuss conditions when such duality is applicable. 

The obtained results can be used in a twofold way: firstly, we show how the infrared limit 
for the chain of oscillators can be described by the corresponding FSG equation with similar 
solutions; secondly, we propose the replacement of the integration scheme of equations with FD 
by the system of coupled Hamiltonian equations. The latter can be considered as a new way of 
analysis of the equations with FD. 



2 Long-range interaction of oscillators 

Consider a one-dimensional system of interacting oscillators described by the Hamiltonian 



M Jo ^ 



\n — m 



COS 



(1) 



where M is mass and Un is displacement from the equilibrium. The last two terms characterize 
an interaction of the chain with the external on-site potential which is defined by the periodic 
function with period a and amplitude J2. The second term in Eq. ^ takes into account the 
interaction of the oscillators in the chain. Here, we consider nonlocal coupling given by the 
power-law function. Constant a is a physical relevant parameter. Some integer values of a 
corresponds to the well-known physical situations: Coulomb potential corresponds to a = 0, 
dipole-dipole interaction corresponds to a = 2, and the limit o; 00 is for the case of nearest- 
neighbor interaction. 

From Hamiltonian it follows the equation of motion. 



1 



m=—oo 



\n — m 



1+a 



Urn + JlUn + J2 Sm {Un) = 0, 



(2) 



where we put M = 1 and the period a = 2tt. 

Let us sketch the derivation of continuous medium equation. We define the field u{k, t) on 
[-K/2,K/2] as 



,t) = ^ Un{t) e 



J^A{Un{t)}, 



where x„ = nAx, and Ax = 2tt/K is distance between oscillators 



1 f 

^«(*) = Kj 



K/2 



dk u{k, t) e 



J^^'{u{k,t)}. 



(3) 



(4) 



These equations are the basis for the Fourier transform which is obtained by transforming 
Unit) = u{nAx,t) from discrete variable Xn = nAx to a continuous one in the limit as Ax — > 
(K oo). The Fourier transform is a generalization of ©, (01) in the limit as Ax 0. Replace 
the discrete Un{t) with continuous u{x,t) while letting Xn = nAx = 27m/K — >• x. Then change 
the sum to an integral, and Eqs. Q become 



u{k, t) 



dx e '"'''u{x,t) = J^{u{x,t)}, 



(5) 
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1 r+co 

u{x, t) = — dk e*^^n(A;, t) = T-^{u{k, t)}, (6) 

where 

u{k,t) = Cu{k,t), u{x,t) = CUnit), Unit) = u{Xn,t), (7) 

and C denotes the passage to the hmit Ax {K oo). The procedure of the replacement 
of a discrete model by the continuous one is defined by the following operations: 

1) The Fourier series transform: 

Ta ■■ Unit) J^A{unit)} = uik,t); (8) 

2) The passage to the limit Ax — > 0: 

£ : u{k, t) C{u{k, t)} = n(A;, t); (9) 

3) The inverse Fourier transform: 

J^-^ : uik, t) J'-^iuik, t)} = uix, t). (10) 

The operation T = L Ta can be called a transform operation (transform map), since it 
performs a tranfrormation of a discrete model of coupled oscillators by the continuous medium 
model. In a similar to jT^ I14j way, we can obtain from @ the equation for n(/c,t) using Q 

d'^uik t) 

' + Jo Mk) uik,t) + Ji uik,t) + J2 ^A{sin(n„(t))} = 0, (11) 

where 

+00 



J — \ ^ ^—iknAx _ \ ^ I —iknAx , AknAx\ _ r„- ( AkAx\ i r „■ /'„- 

Jaik)- 2^ e nilT^ - 2^ (,«^ +e J-Ln+„(e )+M+„(e 



n=— 00 ' ' n=l 



(12) 

and Lii^ai^) is a polylogarithm function. Using the series representation of the polylogarithm 

m 

Li^(e^)=r(l-/3)(-z)'3-i + V^^^^z", \z\<27T, 1,2,3..., (13) 

we obtain 

4(fc) = a, |Axr |fer + 2y;^^i±^-^(Ax)2"(-A:2)", |A:|<1, a ^ 0, 1, 2, 3..., (14) 

where Jq:(0) = 2^(1 + a), C is the Riemann zeta-function and 

a, = 2 r(-a) cos (^) . (15) 

Function Ja{k) can also be presented in the form 

Uk) = 2t^^^. (16) 

n=l 
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from which it follows that 

Ja{k + 2TTm/Ax) = Ja{k), (17) 

where m is an integer. For a = 2, Ja{k) is the Clausen function Cl2{k) |32j . 
Combinig all expressions, Eq. takes the form 

^^^1^ + Jo a„| Axr W uik, t) + 2 Jo ^^^^r^-^iAxf^^i-k'Tuik, t)+ (18) 

n=0 

+ Ji n(A;, t) + J2 .FA{sin (n„(t))} = 0. 
We will be interested in the limit AA; — > 0. Then the Eq. (|18j) can be written in a simplified way 

Q-2u{k, t) + Jo r„,A(A;) u{k, t) + Ji u(A:, t) + J2 J^A{sin (u„(t))} = 0, a / 0, 1, 2, (19) 
where Jo = Jo| Axl™"^"'^} and 



f -|Ax|2-C(a-l)A;2^ 0<a<2, (a / 1) 

"'"^ ^ |Ax|"-VlA:|" -C(a- l)'t^ 2 < a < 4, (a / 3). 



The expression (|2()|) was obtained in jl3| I14j in a slightly different way. Let us note that (|2()j) 
has a crossover scale for 

feo = |aa/C(«-l)r/^'-"^|Axri (21) 

such that Ta^A{k) ~ A;^ for q > 2, A; ^ /co and nontrivial expression a (A;) ~ |A;I" appears only 
for a < 2, A; <C /cq. The crossover was considered also in |141 126j . The expression for Ta^Aik) 
can be considered as a Fourier transform of the operator 

+00 

Ta,nUn{t)= |n-mr(i+") 7x^(t). (22) 

m=— 00 

Performing the transition to the limit k <^ k^ (or more precisely A: Ax <^ koAx), and applying 
inverse Fourier transform to (|19j) gives 



where 



52 

-^u{x,t) + Jo Ta{x) u{x,t) + Ji u{x,t) + J2 sin {u{x,t)) = 0,1,2,..., (23) 



1 - I af let 1 < a < 2, (a 7^ 1) 

r„(x) = T-'{%,{k)} = { ^, _ , ; : J (24) 



C(a-l)^, 2<a<4, (a / 3); 



7; (A:) 



-^1 

'a„|A;|°, 0<a<2, (a / 1) 

-C(a-1)A;2, 2<a<4, (a / 3). 



Here, we have used the connection between the Riesz fractional derivative and its Fourier trans- 
form I2H]: 

aa pp. 



Six 
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The properties of the Riesz derivative can be found in pHl QUI QUI 
If we consider a system with Hamiltonian H = Hi + Hint, where 



H^ 



E 



H 



int 



^ +00 

E 



1 



2 — ' \n — m 

n,m=—oo 



iq:^ 9{un) gium), (26) 



where g{u) is some function of u. Equation Q appears for g{u) = u. The corresponding 
generaUzation of is easily obtained as 



d'^u{x, t) 

at2 



y'(n(x,t)) + Jo 5'(n(x,t))T„(x) (a / 0,1,2, ...), 



(27) 



where g'{u) = dg{u)/du and '7^(x) is the same as in (|^ . 



3 Numerical methods 

We consider + 1 equally spaced by Ax = 2L/{N + 1) coupled oscillators on the finite interval 
{—L,L). The system is described by the equations of motion 

d^Un ^1 , N , ^ 

+ Jo 2^ I _ ^\i+a '^^ + '^1^" + J2smn„ = 0, (n = -y, ...,0, y) (28) 

m=-N/2 ' ' 

and is even. For N ^ 00 Eq. (|28|) coincides with Eq. ((21). Specifing appropriate initial and 
boundary conditions, Eqs. ()28|) could be solved numerically for example by the discretization 
scheme or with the Runge-Kutta method. 

For numerical solutions of the fractional sine-Gordon equation H23|) with < a < 2 

- 

■7-ku{x, t) - Jo Oa— — — ^(x, t) + Ji u{x, t) + J2 sin (ii(x, t)) = 0, (29) 

we have used two methods which are described as the following. 
Method (a): The Riesz derivative could be represented as 07] 

J^u{x, t) = — - {V%u{x, t) + P^n(x, t)) , (30) 

a|x|" 2cos(7ra/2) ^ 

where are Riemann-Liouville left and right fractional derivatives defined by |28| I29| I30| 13 Ij 

1 d"^ r u{i,t)d^ 



V^u{x,t) 



\a—m+l ' 



r(m — a) dx"^ J-oo (^ ~ 

-I)'" d"" Z{i,t )di 
T{m - a) dx"^ (^ - xY 



V-uix,t)= '~'' / , Z^^;'^^"^^, , (31) 



where m — l<a<m. Since we seek for a numerical solution on a finite interval, we will also use 
the Riemann-Liouville left and right fractional derivatives defined on a finite interval (— L, L) 



r(m — a) dx^ J_i (x — ^) 



a—m+l ' 
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Figure 1: Kink-like structures obtained from the evolution of the system of coupled 
oscillators (left panel) and from the solution of the FSG equation (right panel). Bottom 
panel: Comparison of the time evolution of the central oscillator uo{t) (solid curve) with 
the time evolution of the u{0, t) (circles) obtained from the solution of the FSG equation. 
Exponent a = 1.21. 
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Figure 2: Left panel: Breather-like structures obtained from the evolution of the sys- 
tem of coupled oscillators. Right panel: solution of the FSG equation. Bottom figure: 
Comparison of the time evolution of the central oscillator Mo(t) (solid curve) with the 
time evolution of the u{0,t) (circles) obtained from the solution of the FSG equation. 
Exponent a = 1.21. 
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As the next step we approximate the Riesz fractional derivative as follows 

t) ~ {-lV^u{x, t) + a^Vlu{x, t)) . (33) 



' ' ' 2cos(7ra/2) 

To compute Riemann-Liouville fractional derivatives we have used Grunwald-Letnikov dis- 
cretization scheme |28| 

1 

-lT>"u{x, t) = j^^Wk u{x - qh, t), 

q=0 

1 

^Vtu{x,t) = —Y,Wku{x + qh,t), (34) 
g=0 

where h is the discretization parameter, = [{x + L)/h], = [(L — x)/h], [i^] means integer 
part of v, and 

Wk = , . , r. (35) 

r(i + g) r(-a) 

However, for numerical calculations it is more convenient to use, instead of Eq. (|35|) . the following 
recursion relation 

Wk= (^^- ^'^-i' = ^- 

The numerical scheme defined by the Eqs. (jHU could be unstable. Therefore, to obtain a 
stable numerical scheme, one can use shifted Grunwald-Letnikov formulas proposed in j46j 



LV^u{x,t) = l.J2wku{x-{q-l)h,t). (37) 



q=0 



1 

r:Vlu{x, t) = j^^^Wk u{x + {q - l)h, t). 

In our numerical simulations both numerical methods reproduce same results. To ensure the 
stability of the integration scheme the analog to the Courant-Friedrichs-Lewy stability condition 
should be fulfiled 

At/iAx)" < 1/2. (38) 

After substitution of Eqs. (|5H). into Eq. I|29() . we have solved them using the 

explicit finite differences method for the second order time derivative 

52 , , u{x,t + At) + u{x,t- At) -2u{x,t) 

Qf2-(-'i) = • (39) 

We have solved Eqs. (|29|) also by the forth order Runge-Kutta method in time. Results of both 
methods appear to be the same. 
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Method (b): The second method for numerical solution of the FSG equation is based on the 
fact that the Riesz space-fractional derivative admits the explicit representation in the form of 
an integral |2H1 with the limits of integration from zero to infinity. To use this definition for 
the solution of the equation on the finite interval we cut off the upper integration limit and 
approximate the Riesz derivative in the following way 

d°' , , r(l + a) /avrx uix + riA) — 2u{xA) + u{x — -nA) , , 



4 Numerical results 

In this section we compare numerical solutions obtained by two different ways: solution of equa- 
tions of motion for the system of coupled oscillators starting from particular initial conditions 
and solution of FSG with the same initial conditions. 

As the first kind of initial conditions we consider the topological soliton (kink) solution of 
the standard sine-Gordon equation 

^n(O) = 4arctan [Kexp(x)] , (41) 

where for the system of coupled oscillators the variable x defines positions of oscillators, x = nAx. 
The parameter k was fixed to k = 0.001. Three other constants in Eqs. (|28|) and (|29j) were 
Ji = 0.2, Jq = J2 = 0.01. The long-range interaction exponent is a = 1.21. We consider 
-|- 1 = 1001 equally spaced oscillators with the distance between them Ax = 2L/{N + 1), 
L = 500 and boundary conditions Un+2L+i{t) = Un{t) + 2it. For the solution of FSG equation 
the discretization parameter h was choosen to be /\x/2. 

The time evolution of this initial function is shown in Fig. ^ The left panel represents 
the solution of the equations of motion for the discrete system of coupled oscillators, while the 
solution of the FSG equation is shown in the right panel. As it is seen from the figure, solutions 
of both systems are very similar to each other. On the bottom panel we compare the behaviour 
of uo(0 (solid line) and m(0, t) (circles). The fractional order of the space derivative is reflected 
in the time oscillations of the initial soliton-like profile which is the exact stationary solution for 
the standard sine-Gordon equation with a = 2. 

As a second kind of profile we consider two-parametric standing soliton solution of the 
standard sine-Gordon equation 



(0) = 4 arctan 



Kcosh.{x) 



(42) 



where u = 1 and k = 0.05. The constant was choosen such that the ratio wave length to the 
distance between oscillators is small. Three other constants in Eqs. H28|) . (|29() were fixed to 
Jo = 0.01, Ji = J2 = 0.1. The number of oscillators A + 1 = 1001, the distance between 
them Ax = 2L/{N + 1) with L = 500, and the periodic boundary conditions were applied 
Un+2L+i{t) = Un{t)- For the solution of the FSG equation the discretization we choose h = Ax/2. 

Solutions of both equations appear in the form of a breather-like structures. They are ploted 
for different values of a in Figs. 121 121 and HI Solutions for the system of coupled oscillators are 
presented on the left panels of Fig. [3 El and while solutions of the fractional sine-Gordon 
equation Eq. are ploted on the right panels. 
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Figure 5: Profiles of tlie breather- like solutions obtained for the system of coupled oscil- 
lators for t = 100 and a = 1.21, 1.51, 1.91 (solid, dashed and dotted curves respectively) 
in normal, the semi-logarithmic and in the double-logarithmic scales. Straight lines cor- 
respod to /(x) ~ l/x^~^°'. 
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From these figures we conclude that in the infrared hmit solutions of both systems are almost 
identical. Another conclution, that one can draw, is the low sensitivity of the core breather-like 
structure on the exponent < a < 2. In Fig. [S] the snapshots of the breather-like solutions 
obtained for the same initial conditions are plotted for the time instant t = 100. The long-range 
interaction mostly contributs to the behaviour of the tails of the breather-like solutions. The 
exponential spatial decay for short distances is modified by the power-law for large distances 
|35| I36j due to non-local interaction ~ 1/x^^". In Fig. El (middle and low panels) the profiles 
of the breather-like solutions are ploted in semi-logarithmic and in double-logarithmic scales. 
The crossover from the exponential decay of the profile for a short distances to the power-law 
decay for a large distances is explicitly seen. With the increasing value of the parameter a 
the cross-over is shifted for longer time with the agreement to the reduction of the power-law 
interaction to the nearest-neighbor one for a oo. With a decrease of a the exponential part 
of the decay shrinks while the power-law decay is getting broad. 

5 Conclusion 

One-dimensional chain of classical interacting oscillators serve as a model for numerous applica- 
tions in physics, chemistry, biology, etc. Long-range interactions, i.e., with forces proportional 
to l/|a;|^^" are important type of interaction for complex media. An interesting situation arrises 
when, for some reasons, the power a is non-integer. A remarkable feature of these interactions 
is the existence of a transform operator that replaces the set of coupled individual oscillator 
equations into the continuous medium equation with the space derivative of order a, where 
< a < 2, a 7^ 1. Such transformation is an approximation and it appears in the infrared limit. 
This limit helps to consider different models and related phenomena in unified way applying 
different tools of fractional calculus. 

Periodic space-localized oscillations which arise in discrete and continuous nonlinear systems 
have been widely studied in systems with short-range interactions. Here, the system with long- 
range interactions was considered. The method to map the set of equations of motion onto the 
continuos fractional order differential equation is developed in terms of the transform operator. 
It is known that the properties of a system with long-range interactions are very different from 
short-range one. The method offers a new tool for the analysis of different soliton- and breather- 
type solutions in such systems. 

Vice versa, we hope that in some cases the obtained transform operator can be used to 
improve methods for numerical solutions of equations with fractional derivatives. 

In quantum case the application of the transform operator would be highly interesting for 
example to a discrete nonlinear Schrodinger-like equations with the power-law interactions. 
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